#include "Mesh.h"

int main(int argc, char* argv[]){
  EasyMesh M("data/D");
  int n_nodes = M.get_n_nodes();
  std::valarray<Point<double,2>> Nodes = M.get_points();
  for (int i = 0; i < n_nodes ; i++){
    std::cout << i << " : index = " << Nodes[i].get_index() << ", coord= " << std::endl;
    std::cout << Nodes[i].get_coord() << std::endl;
    std::cout << "boundary_mark = " << Nodes[i].get_boundary_mark() << std::endl;
  }
  /*
  for (int i = 0; i < n_nodes ; i++){
    std::cout << i << " : index = " << M.get_point(i).get_index() << ", coord= " << std::endl;
    std::cout << M.get_point(i).get_coord() << std::endl;
    std::cout << "boundary_mark = " << M.get_point(i).get_boundary_mark() << std::endl;
  }
  */
  int n_grids = M.get_n_grids();
  std::valarray<Vector<int,3>> Grids = M.get_grids();
  /*
  for (int i = 0; i < n_grids ; i++){
    std::cout << i << " :" << std::endl;
    std::cout << Grids[i] << std::endl;
  }
  */
  for (int i = 0; i < n_grids ; i++){
    std::cout << i << " :" << std::endl;
    std::cout << M.get_grid(i) << std::endl;
  }
  return 0;
}
